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ABSTRACT 

We present the projected velocity dispersion profile for the remote (d=233kpc) Galactic dwarf 
spheroidal (dSph) galaxy Leo II, based on 171 discrete stellar radial velocities that were obtained from 
medium-resolution spectroscopy using the FLAMES/GIRAFFE spectrograph at the European Southern 
Observatory, Chile. The dispersion profile of those stars with good membership probabilities is essentially 
flat with an amplitude of 6.6±0.7kms _1 over the full radial extent of our data, which probe to the stellar 
boundary of this galaxy. We find no evidence of any significant apparent rotation or velocity asymmetry 
which suggests that tidal effects cannot be invoked to explain Leo IPs properties. From basic mass mod- 
eling, employing Jeans' equation, we derive a mass out to the limiting radius of (2.7 ±0.5) x 10 7 M Q and 
■ a global mass to light ratio of 27-45 in solar units, depending on the adopted total luminosity. A cored 

\ halo profile and a mild amount of tangential velocity anisotropy is found to account well for Leo IPs 

observed kinematics, although we cannot exclude the possibility of a cusped halo with radially varying 
ph ' velocity anisotropy. All in all, this galaxy exhibits dark matter properties which appear to be concordant 

with the other dSph satellites of the Milky Way, namely a halo mass profile which is consistent with a 
I ■ central core and a total mass which is similar to the common mass scale seen in other dSphs. 

2 ■ Subject headings: Galaxies: kinematics — Galaxies: dwarf — Galaxies: stellar content — Galaxies: 

' structure — Galaxies: individual (Leo II) — Local Group 
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Ursa Minor [Palma et al. 2003]), the outer dSph satellites 
might be expected to be relatively unscathed by tidal per- 
turbations. Galactic tides, if acting, would predominantly 
unsettle the outermost regions near the stellar boundary 
of these systems. Moreover, the Draco dSph, though be- 
ing one of the close Galactic satellites, does not show any 
indication of structural perturbations (Odenkirchen et al. 
2001; Segall et al. 2007). Coupled with the negligible 
ratio of ordered stellar rotational velocity (as a manifesta- 
tion of tides; see Section 3 below) to the random motion 
in dSphs, the observations suggest that these galaxies are 
in fact the most dark matter dominated systems found to 
date. It is nevertheless valuable to compare the observed 
properties of dSphs with alternative models, and scenarios 
in which the dSphs are the remnants of tidally disrupted, 
dark-matter-free stellar systems have been investigated by 
a number of authors (e.g. Kroupa 1997; Kuhn & Miller 
1989; Klessen & Zhao 2002). 

Initially, information about the dSphs' dark matter con- 
tent was mainly gleaned from estimates of their central 
velocity dispersion. Fortunately, in recent years wide field 
surveys have provided the opportunity to analyze stellar 
data covering the entire spatial extent of the dSphs and 
to investigate stellar kinematics even beyond their formal 
tidal limits as derived from King profile fits to their stel- 
lar density distributions (usually some tens of arc minutes; 
e.g., Odenkirchen et al. 2001; Kleyna et al. 2003, 2004; 
Tolstoy et al. 2004; Chapman et al. 2005; Mateo 2005; 
Muhoz et al. 2006; Walker et al. 2006; Wilkinson et al. 
2006 and references therein; Koch et al. 2007a). Such 
studies are essential to reliably distinguish between the 
dark matter dominated and tidal disturbance scenarios. 

The mass-density profiles of low-mass galaxies provide 
valuable information about the spatial distribution of dark 
matter on small scales. It has been suggested that halo 
density profiles, in particular the differentiation between 
cored and cusped halos, can provide an insight into the 
actual nature of cold dark matter (CDM; Read & Gilmore 
2005; Colafrancesco et al. 2007; Gilmore et al. 2007). Nu- 
merical CDM simulations with high spatial resolution pre- 
dict dark matter halos that exhibit a central density cusp, 
p(r) cx r~ a with a=l-1.5 (e.g., Navarro, Frenk & White 
1995; hereafter NFW; Diemand et al. 2005). By contrast, 
observations of the kinematics of massive galaxies are fully 
consistent with flat density cores (e.g., Gentile et al. 2004; 
Spekkcns et al. 2005), which also appear to account for 
most of the available data for dSphs (Lokas 2002; Kleyna 
et al. 2003; Read & Gilmore 2005; Goerdt et al. 2006; 
Sanchez-Salcedo et al. 2006; Strigari et al. 2006; Wilkin- 
son et al. 2006; Koch et al. 2007a, Wu 2007). The pres- 
ence of central cores in the dark matter haloes of dSphs has 
important implications for the properties of dark matter 
as it may imply a maximum phase space density for dark 
matter which can be used to discriminate between com- 
peting dark matter candidates (CDM, Warm Dark Mat- 
ter, self- interacting dark matter, etc.). One should bear in 
mind, however, that in the case of dSphs such inferences 
are subject to low-number statistics and uncertainties in 
the models and do not take account of the dynamical evo- 
lution of the dwarf galaxies' progenitor systems (Grebel et 
al. 2003). 

In the light of all the above, the remote dSph Leo II 



provides an interesting test bench to investigate the dark 
matter content of Local Group dSphs. At its present-day 
Galactocentric distance of 233 kpc (Bellazzini et al. 2005) 
it is unlikely to have experienced severe Galactic tides un- 
less it is on a very eccentric orbit that brings it close to 
the Milky Way. Previous kinematic data for this satellite 
were assembled in the study of Vogt et al. (1995; hereafter 
V95) for 31 stars within the core radius at 2.9' (Irwin & 
Hatzidimitriou 1995). From these data, V95 found a M/L 
of 7 and concluded that Leo II must be embedded within a 
massive dark matter halo, although its dark matter prop- 
erties are not extreme compared to those of other dSphs. 
On the other hand, Coleman et al. (2007), using photom- 
etry from the Sloan Digital Sky Survey, estimate a M/L of 
the order of 100, but predicate this value on the assump- 
tion that Leo II is purely dark matter dominated. 

In this Paper, we shall investigate the dynamics of this 
particular galaxy by analyzing data for a large number 
of targets that reach out to the nominal outer optical ra- 
dius at ~ 9' (Irwin & Hatzidimitriou 1995; Coleman et 
al. 2007) thus enabling us to perform a full kinematical 
analysis of Leo II and to derive detailed information on 
its density distribution, its global M/L ratio and thus to 
answer the question of how well-behaved or extreme this 
galaxy is in terms of its dark matter content. 

This paper is organized as follows: our data, their re- 
duction and the derivation of individual velocities are de- 
scribed in §2. We investigate the extent of Leo IPs ap- 
parent rotation in §3 and determine the galaxy's radial 
velocity dispersion profile in §4. These results allow us 
to derive mass and density profiles in §5, while §6 is ded- 
icated to the question of whether Leo IPs density profile 
shows a cusp or a core, as well as the importance of ve- 
locity anisotropy. Finally, we summarize our findings in 
§7. " 

2. OBSERVATIONS AND REDUCTION 

In the framework of the ESO Large Programme 171.B- 
0520(A) (principal investigator: G. F. Gilmore), which 
aims at elucidating the kinematic and chemical character- 
istics of Galactic dSphs, five fields in Leo II were observed 
using the FLAMES multi-object facility at ESO's Very 
Large Telescope (Pasquini et al. 2002). Within the same 
programme, we also analysed the Carina dSph, which is 
described in detail in Koch et al. (2006). In a first paper 
concerning Leo II that made use of these data we presented 
this galaxy's spectroscopic metallicity and age distribu- 
tions, which were derived using the near-infrared calcium 
triplet (CaT) calibration method (Koch et al. 2007b, here- 
after Paper I). For details on our observing strategy, the 
data reduction and calibration techniques, we refer the 
reader to Paper I. In this section we briefly summarize the 
main steps taken. 

2.1. Target selection and acquisition 

The wide field of view of the FLAMES instrument, with 
a diameter of 25', enables one to cover the whole projected 
area of the Leo II dSph with one single pointing, since its 
nominal tidal radius lies at 8.'7 (Irwin & Hatzidimitriou 
1995). However, in order to achieve good sampling out to 
large radii and to yield more information about the stellar 
content near and beyond the tidal radius, we observed five 
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different (overlapping) fields in several instrument config- 
urations, typically offset by 5' with respect to each other 
(see Tables 1 and 2 in Paperl). 

Our targets were drawn from photometry that was 
obtained by the Cambridge Astronomical Survey 
(CASU; Irwin & Lewis 2001) at the 2.5 m Isaac Newton 
Telescope (INT) on La Palma, Spain. From this data set 
we selected red giant candidates with magnitudes ranging 
from the tip of the red giant branch (RGB) at V~18.5mag 
down to 2.5 mag below the RGB tip, reaching as deep as 
21 mag in apparent V^-band magnitude. The total selected 
number of red giant candidates we could thus observe was 
200. Their positions on the sky are depicted in Fig. 1. 

The observations themselves were carried out using the 
identical strategy and instrumental setups as described in 
Paperl, i.e., we used the FLAMES instrument in com- 
bination with the GIRAFFE multifibre spectrograph in 
"low-resolution" -mode centered around the near- infrared 
CaT at 8550A. This set up provides a full wavelength cov- 
erage of 8206-9400A and a resolving power of R ~ 6500 
per resolution element. Although we aimed at exposing 
each configuration for 6hrs to reach a minimum nominal 
signal-to-noise (S/N) ratio of 20 pixel -1 at our spectral res- 
olution of ~ 6500, the major part of the nights was ham- 
pered by bad sky conditions. Consequently, the median 
S/N achieved after processing the spectra is 12 pixel -1 , as 
defined from the variance in continuum bandpasses sur- 
rounding CaT region, from which we will determine our 
velocities in Sect. 2.3.. 

2.2. Data reduction 

Details of the reduction process are given in Paperl. In 
summary, we used version 1.09 of the FLAMES data re- 
duction system, girbldrs, and the associated pipeline ver- 
sion 1.05 (Blecha et al. 2000). After standard bias correc- 
tion and flatfielding, the spectra were extracted by sum- 
ming the pixels along a slit of width I pixel. The final 
rcbinning to the linear wavelength regime was done using 
Th-Ar calibration spectra taken during the daytime. 

Sky subtraction was facilitated by the allocation of 
about 20 fibers per configuration to blank sky. Because 
the pipeline did not perform reliable sky subtraction at 
the time the reduction was performed, custom reduction 
software was written to subtract sky background. First, 
the sky fibers were median-combined to produce an av- 
erage sky for each field. Next, the sky was subtracted 
from each spectrum by modeling the observed spectrum 
as a second order polynomial plus the average sky, and 
minimizing the integral over wavelength of the absolute 
residual. The regions of the spectrum containing the lines 
of the CaT were excluded from the residual calculation 
because they are not expected to conform to our simple 
model of the spectrum. We estimate the accuracy of the 
final sky subtraction to be on the order of ~ 3%, defined as 
the lcr-dispersion of each sky-subtracted spectrum divided 
by the median of the sky spectrum. 

Our data set was then completed by the co-addition 
of the dispersion-corrected and sky-subtracted science 
frames, weighted by the exposures' individual S/N, and 

11 see http://www.ast.cam.ac.uk/~wfcsur 

12 IRAF is distributed by the National Optical Astronomy Observator 
in Astronomy, Inc., under cooperative agreement with the National Sci 



subsequent rectification of the continuum. 

2.3. Radial velocities and Membership estimates 

In order to separate Leo IPs RGB stars from Galactic 
foreground stars (Wyse et al. 2006), we determined the 
individual radial velocities of each target star by means of 
cross-correlation of the three calcium triplet lines against 
synthetic Gaussian template spectra using IRAF'ePh FX- 
COR task (see e.g., Kleyna et al. 2004). The tem- 
plates were synthesized adopting representative equivalent 
widths of the CaT in red giants. The typical median ve- 
locity error achieved in this way is 2.4 km s -1 . Taking into 
account the possibility that the errors returned by FX- 
COR, which are based on the Tonry-Davis R-value (Tonry 
& Davis 1979), may in fact be under- or overestimated 
(e.g., Mateo et al. 1998; Kleyna et al. 2002), we pur- 
sued two different tests of our achieved accuracy. First, 
we determined a multiplicative constant to the formal ve- 
locity uncertainties by deriving radial velocities from the 
first two of the Ca lines at AA8498, 8452 A and separately 
from the third line at 8662 A. The constant was then ob- 
tained by requiring a reduced \ 2 of unity between these 
measurements. This procedure resulted in a re-scaling of 
the FXCOR velocity errors by a factor of 0.6. A second 
estimate of our velocity errors will be given in Sect. 2.4 by 
comparing our individual measurements to the published 
values from V95. A plot of velocity errors versus apparent 
magnitude and S/N is shown in Fig. 2. 

Given LeoII's systemic velocity of ~ 76kms -1 (V95), 
our velocity data set may contain a number of Galactic 
foreground stars along the galaxy's line of sight. However, 
owing to our selection of target stars whose colors and 
luminosities are consistent with their being members of 
Leo IPs red giant branch and due to the much lower Galac- 
tic field star density, LeoII's velocity peak clearly stands 
out against the Galactic contribution (see Fig. 3; and Fig. 4 
in Paperl). From the Besancon synthetic Galactic model 
(Robin et al. 2003) we would expect no more than 4 Galac- 
tic foreground stars within the our color-magnitude selec- 
tion box and with velocities between 40 and 120 km s -1 , 
thus coinciding with LeoII's intrinsic population (see also 
Paperl). 

An initial fit of a Gaussian velocity peak to the Leo II 
data yields a mean radial velocity and velocity dispersion 
of 78.7kms -1 and 7.6kms -1 . Our sample contains 23 
apparent radial velocity non-members, which deviate by 
more than 5 a from this initial fit and thus were rejected 
from further analyses. The mean heliocentric velocity and 
linc-of-sight velocity dispersion were finally determined via 
an iterative error-weighted maximum-likelihood fit assum- 
ing a Gaussian velocity distribution (e.g., Kleyna et al. 
2002; Koch et al. 2007a). From this we find LeoII's mean 
radial velocity to be (79.1 ± 0.6) km s -1 with a global ve- 
locity dispersion of (6.6 ±0.7) km s -1 . These values are in 
good agreement with the measurements of V95, who deter- 
mined a mean velocity of (76.0 ± 1.3) kms -1 and a central 
velocity dispersion of (6.7 ± 1.1) kms -1 based on 31 high- 
resolution spectra of stars within the core radius at 2.'9 
(Irwin & Hatzidimitriou 1995). If one adopts a conserva- 

;, which are operated by the Association of Universities for Research 
ce Foundation. 
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tive 2cr-cut in the velocity distribution in order to define 
membership of Leo II, 158 red giants are included in our 
sample. Relaxing the criterion to a ±3er-cut yields 171 
radial velocity member candidates. Given the low likeli- 
hood of interlopers in our sample, we will apply the 3cr-cut 
throughout the following analyses. 

One noteworthy feature in the velocity histogram 
(Fig. 3, right panel) is the occurrence of an excess of stars 
at ^85kms _1 , whose velocities are approximately 0.9 a 
higher than the sample mean. This peak deviates by 1.8 a 
(taken as yN) from the best-fit curve. We note, how- 
ever, that there is no apparent spatial correlation of the 
targets in this velocity range, nor is there any particular 
separation in color-magnitude space. In order to further 
investigate the possibility that we might be seeing an in- 
dication of a separate cold feature in our data, we ran an 
extensive number of Monte Carlo tests, where we gener- 
ated normally distributed velocity samples centered on the 
observed mean velocity of Leo II assuming the previously 
derived dispersion, and which were additionally varied by 
the measurement errors. It turns out that ~30% of the 
random data sets exhibit peaks in any one bin that devi- 
ate by as much as the observed 1.8 a so that the observed 
excess of stars is not significant above the 1 <T-level. There- 
fore we conclude that the potential second peak in Fig. 3 
is most likely not a real velocity feature, but may rather 
be due to statistical fluctuations. 

2.4. Comparison with other data 

Out of the 31 stars observed by V95, 28 coincide with 
targets in our sample. Since their data are based on high- 
resolution spectroscopy obtained with the HIRES spectro- 
graph at the Keck telescope, it is worthwhile to employ the 
stars in common between the two data sets to assess the 
quality and consistency of our own measurements. Fig. 4 
shows a comparison of the radial velocities of these 28 
stars as measured in this work to those from V95. An 
error-weighted linear least squares fit to these data yields 
the relation v V95 = (0.98 ± 0.10) «This work + (0.5 ± 7.6) 
with a r.m.s. scatter of 4.5kms~ 1 . Although the latter 
value yields an estimate of the expected uncertainty in the 
velocities, it is clear that the errors of both studies and 
the individual S/N ratios have to be considered before one 
employs a simple re-scaling of our values to match the 
quoted r.m.sQ. Hence, we determined scale factors A in 
three bins of S/N (see right panel Fig. 4) such as to yield 
an agreement between our mean velocity errors and the 

scatter in terms of x 2 — J2 1" 2 T " B work ~^ 9 ^ where the 

A "This work + Cr V95 

sum is over the stars in common. As a result, we rescaled 
our final errors by 1.46 (S/N<20), 1.31 (20 <S/N< 40) 
and 1.18 (S/N>40). This leads to the uncertainties finally 
presented in Table 1. 

Recently, Bosler et al. (2007) have performed a low- 
resolution study of the metallicity distribution of Leo II. 
However, their work did not aim at a kinematical analysis 
of this galaxy and yielded a velocity resolution of 54 km s -1 
per resolution element (at the region of the CaT) and a 
median measurement uncertainty of 8.5kms _1 . Hence, 
and since these authors did not publish individual error 
estimates of their measurements, we did not include the 

13 We note that the HIRES data of V95, despite their considerably hi 
contribute a non-negligible fraction to the observed scatter in Fig. 4. 



37 stars that our observations have in common with their 
set of 74 red giants in the comparison in Fig. 4. Nonethe- 
less, an analogous fit to that described above above yields 
v B 07 = (0.93 ± 0.22) v This work + (12.2 ± 17.5) and a 
r.m.s. scatter of lO.Okms -1 . 

In the light of the respective uncertainties associated 
with each of the published data sets, there is no signif- 
icant systematic deviation present in our measurements 
relative to those from V95 and the results agree to within 
the errors. On the other hand, there is an indication of a 
systematic zero point offset of the velocities of Bosler et 
al. (2007). However, this offset is within the uncertainties 
and can be attributed to the lower resolution of the latter 
data set. 

3. TESTS FOR APPARENT ROTATION 

It is often suggested that the apparent rotation of a 
system when viewed in projection can be considered as 
characteristic of tidal perturbation. N-body simulations, 
(e.g., Oh et al. 1995) indicate that the action of the tidal 
field of the Galaxy on a dSph leads to velocity gradients 
in the outer regions of the satellite. When viewed in pro- 
jection, these dynamical effects will lead to a systematic 
change in the mean velocity along the dSphs's major axis, 
thus mimicking rotation (Piatek & Pryor 1995; Oh et al. 
1995; Johnston et al. 1995; Mateo et al. 1998; Read et 
al. 2006a). Hence, the most efficient way of gathering evi- 
dence whether a particular dSph is being affected by tides 
sufficiently strongly that its kinematics are being modified 
is to look for any sign of apparent rotation in its outer re- 
gions. To date, all but one of the dSphs of the Local Group 
show no significant velocity gradients. The sole exception 
is Carina, which exhibits some evidence of a difference in 
the mean velocity at either end of the major axis, for stars 
beyond the nominal tidal radius (Munoz et al. 2006). 

With its present-day Galactocentric distance of 
(233±15)kpc (Bellazzini et al. 2005), Leo II is not ex- 
pected to have been recently affected by tides, and will 
only have been perturbed in the past if its orbit brought 
it into much closer proximity to the Galaxy. This is con- 
sistent with its observed regular structure, as pointed out 
by Coleman et al. (2007). However, since V95 argue that 
Leo II may not have an unusually high dark matter con- 
tent compared to the inner Galactic satellites, it cannot be 
entirely ruled out that it may have been tidally influenced 
during its evolution despite the presence of a (possibly) 
considerable amount of dark matter. 

In order to investigate the question of tides, we imple- 
mented two tests for a significant indication of rotation 
in the Leo II velocities. First of all, in Fig. 5 we plot the 
run of radial velocity as a function of major axis distance. 
In addition, mean radial velocities have been determined 
in radial bins whose widths were chosen so as to main- 
tain a constant number of stars per bin (lower panel of 
Fig. 5). If tides had in fact significantly perturbed the 
outer regions of Leo II, one would expect to see a radial 
gradient in the velocities. In particular, potential tidal 
distortions of the galaxy's outskirts, e.g., reflected in tidal 
tails, could lead to distinct east- west asymmetries in the 
velocity distributions and one would expect excesses of 

gher spectral resolution than our data, are also of low S/N, which will 
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high and low velocity stars (w.r.t. the systemic mean) 
near the outer boundary on opposite sides of the galaxy. 
It is worth noticing that although 13 of our 200 targeted 
stars lie outside of the King-tidal radius of 8.7', only two 
of these lie within ha of the systemic velocity. Recent ev- 
idence suggests that Leo II may in fact be slightly more 
extended, with a nominal tidal radius of 9.2' (Coleman et 
al. 2007), but as these values rely strongly on the concept 
of fitting a particular functional form to surface photom- 
etry, such differences remain insignificant and will not be 
pursued any further in the present work (see Koch et al. 
2007a for a detailed discussion). In order to ascertain that 
we have not missed any potential gradient or high- or low- 
velocity feature in our data due to our selection criteria, 
we repeated the averaging in each radial bin after apply- 
ing several cuts, i.e., at 3 a, 5 a, and 10er of the sample 
mean. For each cut, the data were rebinned to assure a 
constant number of stars per bin. However, as Fig. 5 im- 
plies, there is no significant radial gradient discernible in 
the kinematics of Leo II for any of the velocity cuts. An 
error weighted least squares fit to our velocity data yields 
a slope of (1.2±0.9) kms' 1 kpc -1 . Further, those stars 
with higher and lower velocities than the systemic mean of 
Leo II are equally distributed across the galaxy and there 
is no apparent indication of any localized kinematic excess. 

To investigate the question of the extent to which any 
apparent rotational signal in terms of velocity gradients 
can be excluded, we proceeded a step further and com- 
puted the mean radial velocity difference of stars on either 
side of bisecting lines passing through each individual tar- 
get. In the presence of a distinct apparent galactic rota- 
tion, such a plot should exhibit a clear sinusoidal signal 
(e.g., Walker et al. 2006; Koch et al. 2007a). These ve- 
locity differences are shown in Fig. 6 versus the position 
angle of the respective bisectors. In fact, there is a peak 
with an amplitude of ~2kms _1 seen in our data, which 
is of the same order of magnitude as the radial gradient 
discussed above. Moreover, the maximum of this distribu- 
tion occurs at a position angle of 16.5° ± 2.4°, where the 
error was obtained through bootstrap resampling. In the 
light of the large uncertainty of Irwin & Hatzidimtriou's 
(1995) quoted systemic position angle, i.e., (12° ± 10°), 
we cannot reject a coincidence between the apparent rota- 
tional signal in our velocity data and Leo IPs minor axis. 
To assess the actual reality of such a kinematic gradient, 
we performed 10 Monte Carlo runs, in which we gener- 
ated random samples of velocities at the fixed sky positions 
of our targets. The velocities for these tests were drawn 
from a normal distribution, assuming Leo IPs mean radial 
velocity and velocity dispersion, and additionally allowing 
for a variation by the observed measurement errors. By 
means of these simulations we find that 87% of the ran- 
dom samples exhibit maximum velocity differences larger 
than our observed value (Fig. 6, bottom panel). Hence, 
the apparent rotational signal we detect is only significant 
at a 13% (0.17a) confidence level. 

As Walker et al. (2006) discuss, apparent rotation can 
also originate from relative transverse motion of a dSph 
due to the bulk proper motion of the system. As there are 
no extant observations of LeoII's proper motion, owing 
to its large distance, we cannot approach the question of 
whether the relative motion between the Sun and Leo II 



would be sufficient to account for the full amplitude of our 
observed rotation. All in all, there is no kinematical evi- 
dence of any significant velocity gradients in Leo II, cither 
due to rotational support or produced by Galactic tides. It 
therefore appears that tides most likely have not affected 
Leo II to any significant degree. 

4. RADIAL VARIATION OF THE DISPERSION 

In order to derive the radial velocity dispersion profile 
of Leo II we employed a maximum likelihood technique, 
which is described in more detail in Kleyna et al. (2004) 
and Koch et al. (2007a). In essence, we radially binned 
our data such that the same number of stars per bin was 
maintained. The binsize was defined so as to ensure a 
sufficiently large number of stars per bin, and we chose 
to include no fewer than 12 data points in each bin. A 
Gaussian velocity distribution centered on the single sys- 
temic mean velocity is then assumed for each radial bin. 
This distribution is convolved with the observational errors 
and also with an additional component, which accounts 
for the Galactic foreground contamination. The latter is 
represented by a uniform interloper velocity distribution 
contributing a fraction fi nt to each individual bin (see eq. 
2 in Koch et al. 2007a). As the velocities of our non- 
member stars also appear to be well fit by a power-law 
distribution with an exponent of —0.6, we also considered 
such an interloper distribution in our determinations of 
the dispersion profile. Nevertheless, it turned out that the 
resulting velocity dispersion profiles are indistinguishable 
in practice and fi nt is typically compatible with zero in all 
bins. Finally, the error bounds on the dispersion were cal- 
culated by numerically integrating the total probability of 
the data set and finding the corresponding 68% confidence 
intervals. Fig. 7 displays the radial dispersion profile thus 
obtained. 

Given the results of Sect. 3, we neglected the effect of a 
rotationally-induced inflation of the dispersion in our cal- 
culations. Another critical issue in the proper dynamical 
study of a stellar system is the influence of any significant 
binary population. The presence of any such component 
in kinematical data leads to an additional, artificial in- 
flation of the observed line of sight velocity dispersion. 
Furthermore, the high-velocity extension of a binary dis- 
tribution increases the deviation from a Gaussian so that 
the assumption of a zero binary fraction yields larger er- 
ror estimates, which may invalidate the assessment of any 
radial information derived from the dispersion profiles. 

However, it has been shown via Monte Carlo simula- 
tions in the past that the influence of binary systems in 
dSphs is in fact negligible, as the purely binary induced 
dispersion tends to be small compared to the overall large 
velocity dispersions found in these galaxies (Hargreaves 
et al. 1996; Olszewski et al. 1996). Furthermore, re- 
peat observations of red giant velocities in the Draco dSph 
(Kleyna et al. 2002) and in Fornax (Walker et al. 2006) 
suggest that the impact of binaries on the measured ve- 
locity dispersion is negligible. These kinematical studies 
did not support an overall binary content larger than 40% 
in those dSphs. Also, Kleyna et al. (2002) note that the 
fraction of dynamically significant binaries in the Draco 
sample amounted to less than 5%. Likewise, Koch et al. 
(2007a) showed for the case of the Leo I dSph that the as- 
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sumption of a larger binary fraction does not significantly 
alter the resulting velocity dispersion profile. Finally, there 
seems also to be no evidence that binaries would have sig- 
nificantly affected Leo IPs internal kinematics, as revealed 
by the Monte Carlo simulations of V95. If we interpret 
observed velocity differences as an indication of binarity, 
there is only one object with a difference in velocity be- 
tween our measurement and that of V95 that is in excess 
of 2cr (taken as the measurement uncertainty). If one was 
to interpret this in terms of binarity (not accounting for 
any selection effects in the data) , this would also suggest a 
low binary fraction in Leo 11. The lack of repeat measure- 
ments of the stars in our data set inhibits an assessment 
of potential binaries. Hence, we will proceed by neglecting 
the presence of any binary population in our calculations. 

The resultant velocity dispersion profile (Fig. 7) is es- 
sentially flat out to the last data point and at most subject 
to fluctuations within the measurement errors. The last 
point, shown at a median distance of 0.5 kpc in Fig. 7, 
includes our outermost observed targets at locations of 
1.5 kpc from the galaxy's center. Overall, the profile is 
well described by a constant value of 6.6 kms -1 , which is 
in agreement with Leo IPs global value derived in Sect. 2.3. 
Hence it appears that Leo II does in fact show the same 
behavior as the majority of the dSphs analysed to date, 
in terms of an essentially flat velocity dispersion profile 
(Wilkinson et al. 2006 and references therein). This sug- 
gests that, similar to the other dSphs, this galaxy is also 
a dark matter dominated system out to the largest scales. 
We will turn to the dynamical implications of this finding 
by deriving Leo IPs mass profile and mass to light ratio in 
Section 5. 

4.1. Correlation with metallicities 

For 52 of the targets around Leo IPs velocity peak, 
metallicities were derived in Paper I. From these it was 
inferred that Leo II exhibits a wide range of metallicities, 
covering at least 1.3 dex. The derived age-metallicity re- 
lation indicates that Leo II experienced constant star for- 
mation over an extended period. Overall, this galaxy ap- 
pears to have formed stars from 15Gyr until 2Gyr ago, 
although its dominant stellar population is of intermedi- 
ate age, around 9Gyr (Mighcll & Rich 1996). When split 
into a metal-poor and a metal-rich population by divid- 
ing at the median metallicity, it was shown in Paper I that 
these stellar components do not differ significantly in their 
spatial distributions, i.e., there is no apparent metallicity 
gradient present in Leo II. Moreover, the stellar ages do not 
appear to show spatial trends. It is then worth noticing 
that the star formation histories derived from deep, but 
central, pencil beam surveys such as the HST studies of 
Mighell & Rich (1996) and Dolphin (2002) are also consis- 
tent with our findings in Paper I from a wide-field survey. 
This again argues in favor of the lack of any consider- 
able population gradient between the center and Leo IPs 
outskirts. In dwarf galaxies with extended star forma- 
tion histories, the younger and/or more metal- rich pop- 
ulations generally tend to be more centrally concentrated 
(e.g., Harbeck et al. 2003). Moreover, different stellar pop- 
ulations in dSphs may also have significantly different kine- 
matics, as in the case of Sculptor (Tolstoy et al. 2004) or 
CVn (Ibata et al. 2006). 



In our data, there is no apparent trend of velocity with 
metallicity discernible (top panel of Fig. 8). In the follow- 
ing, we also separated our velocity data according to their 
previously assigned metallicities. The respective radial ve- 
locity histogram is shown in Fig. 8 (left panel). It turns out 
that the metal poor stars have slightly higher velocities on 
average, but the difference between them and the metal 
richer member candidates is marginal. Using the same 
formalism as outlined in Section 2.3, we find the mean ve- 
locity of the metal poor stars to be (79.1 ± 1.6) kms -1 ver- 
sus (76.9 ± 1.0) kms -1 for the metal rich component. The 
overall velocity dispersions for the two sample are indistin- 
guishable and are consistent with that obtained using the 
entire velocity data set. In addition, the data were divided 
into three radial bins (see Fig. 8, right panel). In this case, 
the metal poor (and presumably older; see Paper I) com- 
ponent has a higher velocity dispersion in each bin. How- 
ever, given the low number of stars per bin (of order ten 
at our chosen binning) and the derived measurement un- 
certainties, any difference in the dispersion between both 
metallicity populations is in fact insignificant. Thus we 
conclude that Leo IPs stellar populations are kinematically 
indistinguishable and well-mixed in phase-space, which is 
fully consistent with the lack of any obvious metallicity or 
age gradient (Bellazzini et al. 2005; Paper I). 

5. MASS, LIGHT AND MASS-TO-LIGHT ESTIMATES 

In order to define a mass estimator that can be readily 
implemented we integrated the Jeans equation (Binney & 
Tremaine 1998, eqs. 4-54 ff.; Koch et al. 2007a) under the 
assumption of spherical symmetry and an isotropic veloc- 
ity distribution. Since our data cover the entire face of the 
galaxy, this procedure has the advantage of exploiting the 
full kinematic information that has been gathered in the 
previous sections as opposed to mass estimates that rely 
on the single value of a central velocity dispersion (e.g., 
V95; Coleman et al. 2007 for the case of Leo II). More- 
over, these methods of "core-fitting" make the assumption 
that the mass distribution follows that of the light (Rich- 
stone & Tremaine 1986; V95), which has been shown to 
be invalid in the case of the Draco dSph (Kleyna et al. 
2002). Nevertheless, a limitation of the simple formalism 
employed in this section is the neglect of the degeneracy 
between velocity dispersion and anisotropy (Wilkinson et 
al. 2002), which persists unless large data sets and/or 
higher moments of the velocity distribution are considered 
(e.g., Binney & Mamon 1982; Lokas & Mamon 2003). We 
will account for this effect in detail in Sect. 6. 

The surface brightness profile required for the mass de- 
terminations was adopted from a fit to the data of Irwin & 
Hatzidimitriou (1995). It is well fit by a Plummer profile 
with a characteristic radius of 2.76' ± 0.09' (see Fig. 9), 
where the error was determined via bootstrap resampling 
(e.g., Mackey & Gilmore 2003). This value is in fact very 
similar to Leo IPs core radius from a King (1966) profile 
fit (Irwin & Hatzidimitriou 1995; Coleman et al. 2007) 
and use of the latter value yields analogous results for the 
mass computations. As implied by the results from the 
previous sections the dispersion profile was assumed to be 
radially constant. 

The resultant mass profile of Leo II is shown in Fig. 9. 
From this we estimate its total mass out to the tidal radius 



Stellar kinematics in Leo II 



7 



at 8.7' (0.6 kpc) to be (2.7 ± 0.5) x 10 7 M Q . The quoted 
errors were obtained from a Monte Carlo simulation of the 
mass calculation and take into account the la confidence 
interval of the velocity dispersion as well as the uncertainty 
in the parameters of the light profile. 

Also shown in Fig. 9 is the density profile that we ob- 
tained from the above dynamical considerations. It is 
characterised by a central density of po = (3.4 ± 0.4) x 
10 8 M Q kpc -3 , which is in good agreement with the value 
of (4.0 ± 1.1) x 10 8 M Q kpc- 3 derived by V95 from their 
analysis of LeoII's very center LeoII's density distri- 
bution resembles a cored profile in the innermost regions, 
although we note that our innermost target star lies at 
a distance of 0.5' (0.03 kpc) from the center so that any 
assessment of the "central" density has to rely on an ex- 
trapolation. The profile reaches a power-law slope of — 1 
at approximately 0.13 kpc and converges towards a slope 
of —2 for distances larger than 0.5 kpc. 

With these mass estimates at hand, we derive the 
satellite's M/L ratio by adopting its luminosity from 
the literature. Thence the (M/L) (in solar units) 
varies from 26.6±^ for L v = (9.9±0.2)xl0 5 L Q (V95), 
35. 6±f 8 7 (L y = (7.4±2.0)xl0 5 L Q ; Coleman et al. 2007) 
to 45.41^ for L y = (5.8±1.8)xlO 5 L (Mateo 19980. 
In all cases, a M/L of this order of magnitude strongly 
suggests that Leo II is governed by a considerable amount 
of dark matter on all spatial scales. Mateo et al. (1998) es- 
timate that the stellar (M/L)y for a galaxy like Leo II with 
prominent old and intermediate age populations (see also 
Paper I) is of the order of 1-1.4. If, on the other hand, 
Leo II contained only an old, single-age globular cluster- 
like stellar population, then Mateo et al. (1998) indicate 
that the M/L ratio of Leo II would be higher than the 
value of the galaxy's actual population by a factor of ^1.3. 
Clearly, these results and our observations show that the 
dark matter content of Leo II is well above that expected 
in a purely old stellar population. From their photometric 
structural analysis of Leo II, Coleman et al. (2007) derive 
its M/L via two assumptions: for the case of mass fol- 
lows light they end up with a (M/L)y of 7, whereas for 
the other extreme of a pure dark matter dominance, this 
value may be as high as ~125. The fact that our full kinc- 
matical analysis points towards the higher end of these 
directions underscores the dark matter dominated state of 
Leo II, with at most mild tidal perturbation. Since the 
spatial resolution of our data and the resulting method is 
an improvement over the traditional core-fitting methods, 
which assume a constant M/L, we show in Fig. 9 (bottom 
panel) also the theoretical M/L profile over the full range 
covered by our observations and using the previously de- 
rived parameters. Already at LeoII's center the M/L is 
an order of magnitude higher than the purely stellar esti- 
mates. 

Given the above results, Leo II fits well remarkably well 
on the (M/L) versus My relation for dSphs (see Fig. 5 
in Gilmore et al. 2007), which is an oft-used argument 
in favor of all dSphs being embedded in dark halos of the 
same mass. Although the exact value of such a common 
mass scale is subject to various assumptions such as the 



M/L of the pure stellar component, this underlying halo 
tends to be of the order of 3xlO 7 M (Mateo et al. 1993; 
Wilkinson et al. 2006; Gilmore et al. 2007). 

6. VELOCITY ANISOTROPY - CUSP VERSUS CORE 

All the mass computations in Sect. 5 relied on the 
assumption of an isotropic velocity tensor, i.e., the 
anisotropy parameter [3 — 1 — (vg)/(v,.) was assumed to 
be zero. Nonetheless, there is compelling evidence that 
the shapes and kinematics of most dSphs necessitate non- 
negligible amounts of velocity anisotropy on all spatial 
scales (e.g., Lokas 2002; Kleyna et al. 2002; Wilkinson 
et al. 2004; Klimentowski et al. 2007). As the neglect of a 
non-isotropic velocity distribution can considerably affect 
estimates of the underlying dark halo mass, we explore in 
the following the possible parameter space of Leo IPs veloc- 
ity anisotropy. In doing this we follow the same formalism 
as outlined in Koch et al. (2007a) (see also Lokas 2002). 
That is, under the assumption of a (constant) anisotropy 
parameter and for different plausible mass profiles, we nu- 
merically determined a theoretical velocity dispersion pro- 
file; the model parameters were then optimized in a least- 
squares sense to yield the best representation of the obser- 
vations from Sect. 4. 

6.1. Density cusps 

Motivated by the cosmological simulations of Navarro, 
Frenk & White (1995; NFW), the density profiles of dark 
matter halos are often described by a convenient function 
with only one free parameter, the "characteristic density" . 
For simulated haloes, this parameterisation has proven ap- 
plicable to a wide range of halo masses from the smallest 
scales, such as the dSphs' dark matter associations, up 
to the large scale structures of galaxy clusters. We first 
adopted this cuspy halo mass model, 

M (r) oc M v (ln(l + cr/r v ) - ^ } , (1) 
L 1 + cr/r v J 

in our computations of an empirical dispersion profile. In 
our representation, the only free parameter is the virial 
halo mass, M v , and all other parameters, namely virial 
radius r v and concentration c, are empirically scaled with 
this mass according to the relations obtained by Jing & 
Soto (2000) (see Lokas & Mamon 2001; eqs. 9-12 in Koch 
et al. 2007a for details of the parameterisations). 

The best-fit results for varying degrees of anisotropy are 
overlaid on the observed dispersion profiles in Fig. 10. For 
all choices of (3, the virial mass of the NFW-halo is of the 
order 2.3xl0 8 M Q , corresponding to a formal virial radius 
of about 12 kpc (which is much larger than the stellar lim- 
iting radius, allowing for the factor of ~10 in mass out to 
the virial radius). As for the anisotropy parameter, the 
reduced x 2 f° r a /3 of (—1, —0.5, 0, 0.5) amounts to (0.29, 
0.21, 0.22, 0.47). The best coincidence of prediction and 
observation is obtained when a slight amount of tangential 
anisotropy (/3 ~ —0.25) is included. 

6.2. Cored density profile 



14 Coleman et al. (2007) yield a lower value of po = (1-3 ± 0.4) X 10 s Mq kpc 3 under the assumption of a Gaussian mass and light distribution. 

15 For the subsequent discussions of Leo II's M/L, e.g., Gilmore et al. (2007), we adopt the latter value using Mateo's (1998) luminosity estimate, 
since this value is, in turn, based on the measurements of Irwin & Hatzidimitriou (1995), whose surface brightness profile we employed. 
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The majority of presently available observational data 
of low-luminosity galaxies indicates that their kinemat- 
ics and structures are consistent with cored density pro- 
files (e.g., Lokas 2002; Kleyna et al. 2003; Salucci et al. 
2003; Sanchez-Salcedo et al. 2006). Hence, we also com- 
puted dispersion profiles in the same manner as above for 
a halo profile with variable inner logarithmic slope (Hcrn- 
quist 1990: Read & Gilmore 2005; eq. 13 in Koch et al. 
2007a): 

(7-3)/a 



p(r) cc C 
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As a comparison with extant data of the density profiles 
of Local Group dSphs indicates, the majority of these sys- 
tems are consistent with cored inner halo profiles, which 
is similar to our findings for Leo II (Wilkinson et al. 2006; 
Gilmore et al. 2007 - see their Fig. 4). While the outer 
density profile of Leo II is similar to those of the other LG 
dSphs analysed to date, Leo II apparently exhibits a mildly 
steeper central density profile than the other dwarfs, due 
to the apparently smaller size of its central core. 

From their photometric analysis Coleman et al. (2007) 
argue that Leo IPs M/L ratio has an upper limit of 125 



+56 



under the assumption of a constant density core out to 
their tidal radius. Although our value derived in the pre- 
vious Section is lower by a factor of approximately 3, we 
note that the core size we obtain in the density profile is 
correspondingly smaller. 

7. SUMMARY 

We have obtained a large data set of radial velocity mea- 
surements for about 171 red giant member candidates in 
the remote Leo II dSph, thereby increasing existing pub- 
lished samples more than fivefold (V95). These data ex- 
tend in radial distance out to the galaxy's limiting radius, 
which is an important prerequisite for kinematic studies, 
since Leo II, as the second most remote dSph satellite of 
the Milky Way, is valuable for an assessment of the poten- 
tial role of Galactic tides. With its moderate present-day 
systemic velocity, Leo II can be considered as bound to 
the Galaxy and yet there is no indication of any tidal per- 
turbation. This lack of perturbations is expected if this 
dSph has maintained its present large distance over long 
time-scales. Our study precludes any apparent rotation of 
LeoII's stellar component; moreover we do not detect a ve- 
locity gradient, nor any asymmetry in the sense that high 
and low velocity outliers preferentially reside in particular 
directions along the major axis, which might reflect the 
action of tides. Although other dSphs may well exhibit 
"extratidal" material (e.g., Carina; Muhoz et al. 2006) 
or be in an advanced state of dissolution (e.g., Sagittar- 
ius; Ibata et al. 1994), these effects arc clearly not seen 
in Leo II. Thus we must conclude that this galaxy is a 
purely pressure-supported system which has not been sig- 
nificantly affected by tides over the course of its evolution. 

The radial velocity dispersion profile and resulting den- 
sity profiles show a resemblance to the respective profiles 
of the majority of the other Local Group dSphs, and dis- 
play the well-established features of these systems, namely 
a dispersion which is flat out to the reach of the present 
data and a density profile that is consistent with a cored 
halo mass distribution. We note, however, that our present 
data do not unambiguously distinguish between the pos- 
sibilities of a cusped or cored density profile in terms of 
the x 2 -statistics of the best-fit to our observed dispersion 
profile. Both cored and cusped halo models favor the pres- 
ence of a mild amount of tangential velocity anisotropy in 
the central regions of the galaxy. 

All in all, we find that Leo II is a typical dSph in terms of 
the values of its total mass of ~3xl0 7 M Q and of its char- 
acteristic central density of the order of 3.4 x 10 8 Mq kpc -3 

16 The assumption of velocity isotropy would generally lead to an overestimate of the mass at large radii, if Leo II were indeed tangentially 
anisotropic, whereas if it were radially anisotropic in the outer parts, we would underestimate the mass assuming isotropy. 

17 We note that a radially varying tangential anisotropy /3(r) = r 2 /(r 2 — r 2 ) (Merritt 1985; his model II) yields essentially the same results. 



The parameters were then determined such as to match 
the observed dispersion profile for different values of the 
anisotropy (3. 

As it turns out, the "best-fit" for all choices of 
anisotropy was obtained for a slope parameter of 7 ~ 0, 
which corresponds to a pure core, as opposed to the cuspy 
halo class (7 = 1). Moreover, a scale radius r s of 0.3 kpc 
and a large smoothness parameter a (which determines 
the sharpness of the transition toward the density profile 
at large radii) of ~3.7 were found to yield the best repre- 
sentation of our data (see lower panels of Fig. 10). For this 
cored model, the values of the reduced % 2 are (0.21, 0.20, 
0.21, 0.31) for constant anisotropies of (—1, —0.5, 0, 0.5) 
and, as for the cuspy halos, our data ideally require some 
tangential anisotropy of the same order of magnitude, that 
is 0.25. 

It is worth noticing that our data do not allow us to un- 
equivocally constrain the shape of the velocity anisotropy 
tensor at large radii for either choice of halo profile. How- 
ever, we can conclude from Fig. 10 that strongly radial 
anisotropy is not suitable to account for the dispersion 
profile towards the central regions. 

6.3. Cusp or Core ? 

For both profiles investigated above, the mass at the 
galaxy's limiting radius amounts to ~ 2.0 x 10 7 M Q and 
is thus slightly smaller than that derived from a purely 
isotropic velocity distribution in Sect. e Pi . 

Judging by the values of the x 2 statistics, a cored den- 
sity profile is favored to account for the Leo II observed 
kinematics; however, the difference between this and the 
cuspy model is small. In particular, there is an indication 
of a leveling off of the predictions for the dispersion to- 
wards larger radii for the cored model, a trend that merits 
further investigation from a large sample of stars at and 
beyond the limiting radius of this galaxy. 

In addition, we re-ran the computations with a ve- 
locity anisotropy which varies radially according to the 
prescription of Osipkov (1979) and Merritt (1985), i.e., 
f3(r) — r 2 /(r 2 + r 2 ), which yields radial anisotropy for 
the case of r > r a . The resulting curves are shown in 
the right panels of Fig. 10. While the NFW halo requires 
an anisotropy radius r a of ^0.7 kpc and thus larger than 
LeoII's tidal radius, the cored profile is best associated 
with a r a of ~0.2kpc. The respective \ 2 statistics are 
then 0.22 and 0.19, thus again marginally favoring a cored 
haldHI 
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(Gilmore et al. 2007). Depending on the adopted to- 
tal luminosity of Leo II, we estimate its M/L ratio to lie 
in the range 25 to 50. This is, in conjunction with the 
flatness of the dispersion profile, striking evidence that 
Leo II is a dark matter dominated system. Moreover, the 
value we find is consistent with the idea that all dSphs are 
embedded in a dark matter halo with one single charac- 
teristic mass scale. Taking into account that all nearby 
dSph galaxies contain very old stellar populations (Grebel 
& Gallagher 2004) , it thus appears reasonable to consider 



these objects as the smallest surviving condensations of 
dark matter with associated stellar populations in the uni- 
verse (Read et al. 2006b). 
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Fig. 1. — Distribution of our targets on the sky, centered on Leo II. Apparent radial velocity non-members are designated with open circles. 
The solid and dashed lines indicate LcoII's core and nominal King tidal radius at 2.'9 and 8f7 (Irwin & Hatzidimitriou 1995). 
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Fig. 2. — Velocity errors versus V-band magnitude (left panel) and Signal to Noise ratio per pixel (right panel). The dashed lines in 
the right panel indicate our formal mean velocity errors in each S/N bin, whereas the solid lines show the r.m.s. discrepancy between our 
measurements and those of V95 (open squares), which we use to re-scale our uncertainties in Sect. 2.4. 
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Fig. 3. — Histograms of the radial velocities of stars in our survey of Leo II. The left panel displays the full sample, whereas the right panel 
only shows stars around Leo IPs systemic velocity. The best-fit Gaussian to this distribution is shown as a solid line. 
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Fig. 4. — Comparison of our velocity measurements to those estimated by Vogt ct al. (1995) for the 28 stars in common between the two 
data sets. The scale factors to the error bars have not been applied in this diagram to illustrate our internal formal measurement errors. The 
dashed line is unity. 
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Fig. 5. — Top panels: Radial velocity versus major axis radial distance. The right panel shows those stars around Leo IPs mean systemic 
velocity, while the full data arc included in the left panel. Members within (outside) a 3<r cut are shown as filled (open) symbols. The bottom 
panel displays mean radial velocities obtained in radial bins with a constant number of stars, where we adopted different cuts in velocity to 
consider the presence of high- or low- velocity member stars: 3a (open circles), 5cr (solid squares), and 10 a (filled diamonds). No significant 
radial gradient is discernible in our data. Monte Carlo tests show that the hint of a change in the mean velocity between the outer bins at 
each end of the major axis is not statistically significant (see text for details). Note the different scales in each of the subplots. 
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Fig. 6. — Test for apparent rotation in our observations (top panel) and a sample of Monte Carlo simulations (bottom panel). The points 
with error bars show the mean velocity difference of stars on either side of bisecting lines at the respective position angles. The best-fit 
sinusoid is displayed as a solid line. The discernible rotation amplitude at a position angle of ~23° is statistically insignificant — 87% of the 
random samples as shown in the bottom panel produce a comparable amplitude or larger. Indicated as a dashed line in the top panel is the 
minor axis position angle of 12° from Irwin & Hatzidimitriou (1995). 
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Fig. 7. — Leo IPs radial velocity dispersion profile obtained for a 3 a membership criterion. In each panel, each bin contains the indicated 
constant number of stars. 
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Fig. 8. — Upper panel: Velocity distribution for a subsamplc of our data for which there is extant metallicity information (Paper I). The 
velocity histograms and dispersion profiles in the lower panels have been separated into a metal-poor and a metal-rich component, which 
exhibit only marginally different kinematics. 
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Fig. 9. — Mass (middle left panel) and density (middle right) profile of Leo II, obtained from Jeans equation under the assumption of 
an isotropic velocity distribution and a constant velocity dispersion. Error bounds originate from the uncertainties in the light and velocity 
dispersion profile. The surface brightness profile entering the calculations is shown in the top panel (after Irwin & Hatzidimitriou f995). The 
solid line is a best-fit Plummer profile. At Leo IPs distance, 1' corresponds to about 0.07 kpc. Based on mass and light profiles, we derive the 
M/L distribution shown in the lower panel. 
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Fig. 10. — Dispersion profiles derived under the assumption of varying degrees of velocity anisotropy j3. The upper and lower panels adopt 
two different density models, namely a NFW cusp (upper panels) and a cored halo-profile (lower panels), respectively. While the left panels 
maintain constant values of /3, the right panels were determined for a velocity anisotropy which varies radially according to the Osipkov-Merritt 
(O.M.) prescription with an anisotropy radius as stated. 



Stellar kinematics in Leo II 



17 



Table 1 

Measured properties of radial velocity members in Leo II 



Star 


a (J2000) 


6 (J2000) 


r['j 


Vr 


°v r 


[Fe/Hj 


a\Fe/H\ 


T_18 


11 13 03.5 


22 11 35.0 


7.06 


48.49 


3.18 


-1.34 


0.18 


TU9 


11 13 23.0 


22 09 23.0 


1.64 


78.18 


5.04 


-1.62 


0.52 


V.19 


11 13 35.6 


22 09 06.0 


1.68 


76.33 


8.59 






T_20 


11 13 32.1 


22 09 12.0 


0.77 


90.10 


16.29 






T_21 


11 13 22.7 


22 08 09.0 


2.11 


83.99 


3.53 


-1.75 


0.25 



Note. — This Table is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance 
regarding its form and content. Metallicities are those from Paper I. 



